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Matter with chiral fermions is microscopically described by theory with quantum anomaly and 
macroscopically described (at low energy) by anomalous hydrodynamics. For such systems in the 
presence of external magnetic field and chirality imbalance, a charge current is generated along the 
magnetic field direction — a phenomenon known as the Chiral Magnetic Effect (CME). The quark- 
gluon plasma created in relativistic heavy ion collisions provides an (approximate) example, for 
which the CME predicts a charge separation perpendicular to the collisional reaction plane. Charge 
correlation measurements designed for the search of such signal have been done at RHIC and the 
LHC for which the interpretations, however, remain unclear due to contamination by background 
effects that are collective flow driven, theoretically poorly constrained, and experimentally hard to 
separate. Using anomalous (and viscous) hydrodynamic simulations, we make a first attempt at 
quantifying contributions to observed charge correlations from both CME and background effects 
in one and same framework. The implications for the search of CME are discussed. 
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Introduction .— The study of matter with chiral 
fermions has generated significant interest recently, en¬ 
compassing a wide range of systems from condensed mat¬ 
ter materials to hot dense nuclear matter [1]. Of particu¬ 
lar interest, are possible anomalous effects that can man¬ 
ifest the microscopic quantum anomaly of chiral fermions 
in the macroscopic transport properties of matter. The 
universal nature of chiral anomaly often leads to certain 
universal features of such anomalous transport effects. A 
well-known example is the Chiral Magnetic Effect (CME) 
the generation of a vector current J (a parity-odd 
vector quantity) along an external magnetic field B (a 
parity-even axial-vector quantity): 

J = CaUaB ( 1 ) 

where /a a is a nonzero axial chemical potential that quan¬ 
tifies the amount of chirality imbalance i.e. the difference 
in numbers of right-handed and left-handed fermions. 
The coefficient Ca is the universal constant originated 
from the chiral anomaly coefficient, e.g. Ca = N c e/(2ir 2 ) 
for each flavor of massless quarks in QCD. 

One concrete physical system where the CME may 
occur and get experimentally observed, is the quark- 
gluon plasma (QGP) — an extremely hot, deconhned 
form of nuclear matter that has been created and mea¬ 
sured in high energy nuclear collisions at the Relativistic 
Heavy Ion Collider (RHIC) and the Large Hadron Col¬ 
lider (LHC) [2-4], Evidently three elements are needed 
for (1) to happen. First a chiral QGP with (approxi¬ 
mately) massless quarks is necessary for anomaly effect. 
While the spontaneous breaking of (approximate) chiral 
symmetry is a fundamental property of QCD vacuum, it 
is indeed predicted by Lattice QCD simulations as well as 


theoretical models that such symmetry is restored at the 
high temperature accessible in heavy ion collisions. Fur¬ 
thermore a chirality imbalance /r a 7^ 0 is needed. This 
pertains to a salient feature of QCD as a non-Abelian 
gauge theory: the topologically nontrivial gluonic con¬ 
figurations such as instantons and sphalerons that are 
known to be crucial for understanding nonperturbative 
dynamics of QCD. These configurations couple to quarks 
through chiral anomaly and “translate” topological fluc¬ 
tuations into chirality imbalance for quarks, thus creat¬ 
ing nonzero /.la on an event-by-event basis. Finally in 
a heavy ion collision, very strong magnetic field results 
from the incoming nuclei that are positively charged and 
move at nearly the speed of light. Such B field has a 
magnitude on the order of eB ~ m 2 and points ap¬ 
proximately in the out-of-plane direction [5-7]. A un¬ 
ambiguous observation of CME in heavy ion collisions 
would therefore provide experimental evidence for a chi¬ 
ral symmetric QGP as well as the QCD topological con¬ 
figurations. In addition to the CME, various other inter¬ 
esting anomalous transport effects have been proposed, 
such as Chiral Separation Effect [8, 9], Chiral Electric 
Separation Effect [10, 11], Chiral Magnetic Wave [12— 
14], Chiral Vortical Wave [15]. For recent reviews see 
e.g. [1, 16-18], 

In this study we focus on the Chiral Magnetic Effect 
in heavy ion collisions. The CME (1) predicts a charge 
separation along the out-of-plane direction with excessive 
positive charges accumulating on one tip of the fireball 
and negative charges on the other tip. Such a separation 
can be measured by the following reaction-plane depen- 
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dent azimuthal correlation observable: 

Ob/3 —O COS ((j>i ~b (fij ^'/(’Rp) ba/3 (2) 

with a, /? = i labeling the charged hadron species 
and 4>ij the azimuthal angles of two final state charged 
hadrons. The >Frp denotes reaction plane angle and for 
later convenience we set d^Rp = 0. This observable has 
been measured at RHIC [19-22] for a variety of beam 
energy and centrality as well as at the LHC [23]. The 
measurements show highly nontrivial change-dependent 
azimuthal correlations, i.e. charge asymmetry is signifi¬ 
cant in high energy collisions and disappears at low en¬ 
ergy. While some aspects of data are consistent with 
CME expectations, an unambiguous extraction of CME 
signal has been obscured due to signifiant background ef¬ 
fects driven by bulk flow [24-29]. This has been clearly 
revealed by examining another correlation observable: 

&af3 —b COs(</>; (f)j ) bct/3 (3) 

for which data show opposite trends from CME expec¬ 
tations. It was found that the transverse momentum 
conservation and the local charge conservation can make 
strong contributions to these observables. For detailed 
discussions see e.g. [17, 18]. 

Given the importance of CME and given the presently 
unclear situation in experimental search, what is criti¬ 
cally needed is a quantitative modeling of both the CME 
signal and the pertinent background effect that would al¬ 
low a meaningful comparison with data. Let us identify 
a number of outstanding challenges faced in such effort: 
(1) a description of CME in the hydrodynamic framework 
that is built on top of state-of-the art, data-validated bulk 
evolution for heavy ion collisions; (2) a quantification of 
the influence of key theoretical uncertainties like initial 
axial charge fluctuations and magnetic field lifetime on 
the CME signal; (3) an evaluation of background con¬ 
tribution consistently in the same bulk evolution frame¬ 
work; (4) predictions for further measurements that can 
help verify theoretical assumptions in the modeling. It is 
the purpose of this Letter to report a significant step for¬ 
ward in addressing these outstanding questions and thus 
substantially advancing the search of CME in heavy ion 
collisions 

CME signal from anomalous hydrodynamics .— The 
Chiral Magnetic Effect (1) implies anomaly-induced con¬ 
tributions to hydrodynamic currents, and a first step one 
needs to take is to integrate CME contribution with the 
usual viscous hydrodynamical simulation of heavy ion 
collisions. The theoretical foundation for this integration 
has been laid down recently. Fluid dynamical equations 
with chiral anomaly, i.e. anomalous hydrodynamics, have 
been derived [30]. ( For out-of-equilibrium situation, see 
[31] in which anomaly effects are incorporated in the 
framework of kinetic theory.) Explorative attempts were 
recently made to apply them for phenomenological mod¬ 
elings in heavy ion collisions [32-34]. In this work, we 


adopt the approach similar to that in [33], which treats 
the fermion currents as perturbations and solves anoma¬ 
lous hydro equations for these currents on top of the 
data-validated viscous hydrodynamic background. The 
equations read: 


= d M K + QfCAVAB^) = 0 (4) 

= d „ ( n A uP + QfCAiivB ") = -Q 2 f eC A E^( 5 ) 


where E ^, B M are covariant form of electromagnetic 
fields. Note these equations are for each quark flavor with 
corresponding charge Qf. The flow field iC and local 
temperature are taken from background hydro solution 
by “VISH2+1”, a 2 + 1 viscous hydrodynamics code as¬ 
suming boost invariance [35]. The quark susceptibility at 
given (local) temperature that relates density with chem¬ 
ical potential is taken from lattice results [36]. The evo¬ 
lution is followed by a slightly generalized Cooper-Frye 
freeze-out procedure (see [33] for technical details) that 
accounts for nontrivial charge transport. Starting from 
an initial condition of nonzero axial charge density, these 
equations indeed lead to a spatial separation of positive 
and negative charges on the freeze-out surface along B 
direction. Combined with strong radial flow this leads to 
an event-wise azimuthal distribution of charged hadrons 
of the following form: 


dN H 

d(j> 


oc [1 + 2Q H di sin (<f>) + ...] 


( 6 ) 


J CME 


where “H” labels the species of the hadron, e.g. H = 
7T^ , AT ,.... The a E , computed from the anomalous hy¬ 
dro equations, quantifies a CME-induced out-of-plane 
charge separation. This gives a contribution to observ¬ 
able (2) and observable (3) by E = -S^ E = 

—Q H °‘Q Hf> {ai a a^ f> ). The quantitative results depend 
on two key input factors in the simulation, which we dis¬ 
cuss below. 

The first is the initial axial charge density that could 
be generated from gluon topological fluctuations which 
is a most significant theoretical uncertainty. A plausible 
strategy is to study its influence on CME signal and to 
put constraint on such uncertainty through data com¬ 
parison. A reasonable assumption is to have initial axial 
density per each flavor of light quarks to be proportional 
to initial entropy density sj, with a proportionality con¬ 
stant Aa = {ua/N f) / si s=s QA/{NfS ) where Qa is the 
total initial axial charge while S the total entropy. We 
note in the linearized regime (owing to the fact that these 
density fluctuations are all small), this smooth average 
initial condition is essentially equivalent to event-wise lo¬ 
calized axial density “lump” with a probability distri¬ 
bution proportional to initial entropy density, and thus 
is not so much different from the event-by-event simula¬ 
tions in [34]. In our modeling the CME signal is found 
to linearly depend on parameter A a, be. a E = A . 
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The second is the magnetic field B. While its peak 
magnitude at initial impact of collision has been deter¬ 
mined [5-7], its subsequent time evolution is affected by 
the created partonic medium and not fully understood 
at the moment [37]. The CME results crucially depend 
upon the lifetime of B field and it is vital to under¬ 
stand such dependence. We take B to be homogeneous 
in transverse plane and use a parametrization eB(r) = 
(eB) 0 /[l + ( t/t b ) 2 ], with (centrality-dependent) peak 
value {eB) 0 from [7], and study how the CME signal 
depends on the lifetime tb- 



FIG. 1. (Color online). (Left) Dependence of af + on tb at 
50 — 60% centrality(6 = llfm). (Right) Dependence of af + 
on centrality for tb = lfm/c and 3fm/c. 


In short the CME signal is controlled by the two key 
parameters: that characterizes initial axial charge as 

well as the magnetic lifetime tb- In Fig.l we show our 
results for such dependence. Note throughout the paper 
we focus on RHIC AuAu collisions at yfs = 200GeV. 

Background contribution from transverse momentum 
conservation. — We now turn to discuss background ef¬ 
fects which have to be understood before any compari¬ 
son with data. Past studies have already revealed that 
the opposite-charge pair correlations are likely dominated 
by the local charge conservation effect that overshadows 
CME contribution [25, 28, 29]. The same-charge pair 
correlations suffer mainly from the transverse momentum 
conservation (TMC) effect which is likely comparable to 
CME signal [25, 27, 38]. Furthermore the way to quantify 
TMC contribution is theoretically well understood [27]. 
For this reason, we will focus only on the same-charge 
pair correlations in the present study. 

Intuitively the TMC-induced correlations can be un¬ 
derstood as follows. Imagine N particles each with 
the same momentum magnitude P along the same axis, 
then any particle’s momentum is to be balanced by all 
the other particles with each balancing on average the 
amount of —P/{N — 1). For N » 1 this implies a 
back-to-back correlation on the order ~ \/N in each 
momentum direction. Therefore the TMC gives nega¬ 
tive contribution ~ — 1 /N to both in-plane projection 
(cos <f>i cos 4> 2 ) and out-of-plane projection (sin (j>i sin <f> 2 ). 
Due to elliptic flow however the effect is slightly stronger 
for the in-plane, and as a result the TMC contributes 
7 ™ c ^ -~y 2 /N to observable (2) and S TMC ~ — 1/N to 
observable (3). 

To quantitatively evaluate this background effect, we 
have first generalized the analytic formula for single¬ 


component TMC in [27] to the case of multiple types 
of hadrons. From that we can compute the TMC con¬ 
tributions to two-particle azimuthal correlations, assum¬ 
ing a single particle distribution of the form Jh(p±) = 
/o,ff(p±)[l + 2 u 2 l jj(p±)cos( 20 )]: 


cTMC 1 TMC . _ [(PJ-)a(l ± V2,a)][(PJ.)/3(l ± V 2 ,p)] 

a0 7a/? N tmc (pI)( 1 ±v 2 ) 

( 7 ) 

Here, by definition of (2), (3), S + 7 and 5 — 7 equal to 
2(cos0cos<^>) and 2(sin<^>sin0) respectively. In the above 
quantities with hadron labels a, /3 are computed for the 
specific hadron species while those without are computed 
from all hadrons. We have also introduced the following 
notations: v 2 ,h = ( P±v 2 ,h)/(p ±), i>2,h = (p 2 ± v 2}H ) / (p 2 ±), 
with (...) = [f^’> dp ± {.. .)p ± f 0 , H ]/[Jp^ dp±p ± f 0 ,H]- 
Here kinematic cuts are ( p_l,<,p_l,> ) = (0.15,1.85) GeV 
in accord with pertinent experimental analysis. The 
A^tmc should be the total number of all produced parti¬ 
cles which is closely related to observed multiplicity but 
not exactly that due to kinematic constraints and detec¬ 
tor efficiency of experiments. Given such practical un¬ 
certainty, we treat A%mc as a parameter that controls 
the magnitude of TMC effect. With the above formulae 
we then compute the TMC contributions to both observ¬ 
ables (2) (3) using the final hadrons’ spectra from the very 
same viscous-hydro calculations as also used for comput¬ 
ing CME signals. 

Understanding data for same-charge pair 
correlations. — Up till now we have built a hydro- 
based framework for simultaneously quantifying both 
CME and TMC contributions to observables, with three 
undetermined model parameters (A ,4 and tb for CME 
while A^tmc for TMC). This allows one to constrain such 
parameters by comparison with data. Note the data 
we compare to, are for all charged hadrons, mainly 7 ^. 
As aforementioned, we focus on the same-charge pair 
correlations for which the following working assumptions 
are very plausible [25, 27, 38, 39]: 


^,data 

7a, p 


-CME 

7a, P 


TMC 
7a,/9 


rdata ?CME 
°a,f3 — °a,p 


I xTMC 
+ °a,/3 


.( 8 ) 


where the two types of hadrons H a ,Hp are either both 
positively or negatively charged. It shall be noted that 
the so-called “dipolar flow” (see e.g. data from STAR [40] 
and ATLAS [41]) may in principle make a contribu¬ 
tion to the above correlations. This effect in the inte¬ 
grated [0, 2]GeV/c region (where STAR measurements 
are made) is negligible, but may be significant in the 
integrated [2,5]GeV/c region. This latter point may be 
important for understanding the ALICE data [23] for 7 , <5 
correlations, measured in the full [0,5]GeV/c region. 

Now by noting 7 CME = —JCME, one 0 f,t a i ns f rom (8) 
that (7 + < 5 )™ c « (7 4 - 5 ) data from which one can fix 
the parameter IVtmc- This fitting gives IVtmc/A part = 
(14,14,15,17, 21) for centrality class (10—20, 20—30, 30— 
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40,40 — 50,50 — 60)% respectively (with 7V part the par¬ 
ticipant number for given centrality), which appear very 
reasonable [27, 38]. With the TMC contribution fixed, we 
can then subtract it from data and determine the magni¬ 
tude of CME signal from (sin (p sin <j>) = (8— 7 )/2 channel: 
see Fig.2 (left). This then allows extraction of initial con¬ 
dition parameter A a after assuming reasonable value for 
magnetic field lifetime tb- In Fig .2 (right) we show the 
extracted A a versus centrality for tb = 1 and 3fm/c re¬ 
spectively. In these plots we have included error bars that 
originate from the uncertainty in the STAR data, which 
is dominated by systematic uncertainty in flow measure¬ 
ments and is on the average level of 15% [19, 20]. 
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FIG. 2. (Color online). (Left) (<5++ — 7++)/2 ( i.e. 
(sin rf>i sin 1 / 2 )++) data and its decomposition into CME and 
TMC contributions for different centrality. (Right) Depen¬ 
dence of extracted initial axial charge parameter Aa on cen¬ 
trality for tb = lfm/c and 3fm/c. 

A few remarks are in order. First, we point out that 
the above scenario is able to provide a description of data 
with modest values of the lifetime of magnetic field. Sec¬ 
ond, CME contribution increases sensitively with larger 
initial axial charge Aa, and even if the actual B lasts 
say just a fraction of fm/c, it may be easily compensated 
by a mild increase of A a- Finally we notice the above 
extracted initial axial charge and its centrality trend 
are consistent with theoretical estimates [42] based on 
Chern-Simons diffusion rate for gluonic topological fluc¬ 
tuations. The estimate in [42] gives dQ^/dp ~ (20 — 40) 
and dS/drj ~ 3600 which gives Aa ~ (0.05 — 0.1 )/Nf, in 
consistency with our results. Furthermore the estimate 
in [42] suggests d(Q\)/dr] ex ( dS/dr ?) 4 / 3 therefore one in¬ 
fers A^ oc ( Q\)/S 2 oc S~ 2 / 3 which implies Aa increases 
mildly from central toward peripheral collisions, indeed 
also in line with our results. 

Predictions for identified particles’ azimuthal 
correlations. — Clearly what have been achieved so 
far, are ( 1 ) the use of a hydro-based framework to estab¬ 
lish quantitative relations between the CME and TMC 
contributions to observables and the model parameters, 
and ( 2 ) the determination of preferred model parameters 
that describe data. As a natural and mandatory next 
step, one needs to go beyond understanding just existing 
data and make unambiguous predictions that can be 
tested with future data. Such test is crucial for verifying 
the present physical interpretation ( 8 ) of same-charge 
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FIG. 3. (Color online). Predicted correlations versus 

centrality for identified hadron channels 7r + 7r + , K + K + , pp, 
AA, and pA for symmetric two-flavor (upper) and three-flavor 
(lower) cases. The computation uses tb = lfm/c. To improve 
readability, the error bar, which is on the average level of 15% 
(see text and Fig. 2), is not included. 

We propose to use the same-charge azimuthal correla¬ 
tions for various identified hadron species as a nontrivial 
further test. As evident from (4) the CME occurs for 
each light quark flavor and leads to flavor-specific sepa¬ 
ration in the QGP. Upon hadronization such quark-level 
charge separations convert into hadrons according to 
their quark contents. On chemical freeze-out surface, one 
constructs hadron chemical potentials out of quarks e.g. 

AA+(uci) AC AO, h’K+tus) AA AG, l^p^uud) l^d: 

Afofuds) = AA + AG + AG, etc, where the charge sepa¬ 
ration effects are encoded in the spatial dependence of 
the quarks’ chemical potentials. With our model pa¬ 
rameters already fixed above, we can make quantitative 
predictions in the same framework for various identified 
hadrons’ same-charge correlations (2). In addition this 
allows a possible test of how “chiral” the strange quark 
may be in the QGP, as the effects for strange hadrons de¬ 
pend on whether s quarks have a CME-induced separa¬ 
tion [43]. We have done computations for both two-flavor 
and three-flavor cases and the results for various channels 
are shown in Fig.3, to be tested by future measurements. 

Conclusion. — Anomalous transport effects like the 
CME provide a newly emergent route of probing a chi¬ 
ral symmetric quark-gluon plasma in heavy ion collisions. 
Such global chiral effects rely upon manifestation of chi¬ 
ral symmetry and axial anomaly in macroscopic, hydro- 
dynamic quantities such as vector and axial currents. An 
unambiguous experimental observation of CME is of fun¬ 
damental interest and could be a signature toward iden- 
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tification for the “on” and “off” of QCD chiral restora¬ 
tion in the RHIC Beam Energy Scan experiments. To 
achieve this goal, however, requires quantitative model¬ 
ing of CME and background effects that can be mean¬ 
ingfully validated by data. 

In this study we have used the anomalous hydro¬ 
dynamics framework, combined with usual viscous- 
hydrodynamic simulations of bulk evolutions for heavy 
ion collisions, to quantitatively evaluate the charge sep¬ 
aration induced by Chiral Magnetic Effect as well as the 
background contributions from Transverse Momentum 
Conservation. We have identified the key parameters 
and studied how the resulting experimental observables 
depend upon these parameters. A main finding is that, 
with very plausible choices of model parameters, in partic¬ 
ular the magnetic field lifetime and the initial condition 
for axial charge, a successful description of the present 
same-charge pair correlation data can be achieved. Given 
the modeling framework with fixed parameters, nontriv¬ 
ial patterns for identified hadron correlations in various 
channels have also been predicted. This would allow fu¬ 
ture verification of the proposed interpretation of present 
data in terms of CME plus background. As a final re¬ 
mark, the CME contribution to the opposite charge cor¬ 
relations has been computed in our framework, but to 
compare with data would require a very careful quan¬ 
tification of the overwhelming background contribution, 
which will be studied in a future work. 
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